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ABSTRACT 

The total area coverage by starspots is of interest for a variety of reasons, but direct 
techniques only provide estimates of this important quantity. Sunspot areas exhibit 
a lognormal size distribution irrespective of the phase of the activity cycle, implying 
that most sunspots are small. Here we explore the consequences if starspot areas were 
similarly distributed. The solar data allow for an increase in the fraction of larger 
sunspots with increasing activity. Taking this difference between the size distribution 
at sunspot maximum and minimum, we extrapolate to higher activity levels, assum- 
ing different dependencies of the parameters of the lognormal distribution on total 
spot coverage. We find that even for very heavily spotted (hypothetical) stars a large 
fraction of the spots are smaller than the current resolution limit of Doppler images 
and might hence be missed on traditional Doppler maps. 
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1 INTRODUCTION 

In recent years, an ever increasing number of spotted stars 
have been mapped using Doppler imaging. The maps reveal 
the surface distribution of starspots, which in general are 
large compared to even the largest sunspots. While Doppler 
images do a good job of catching starspots that modulate the 
line profile, it is extremely difficult to detect a background 
of small starspots more or less homogeneously distributed 
over the stellar surface. While TiO-band mapping still suf- 
fers from notable uncertainties (see Sec. 2), it should in prin- 
ciple be able to pick up non-modulating and homogeneous 
spot distributions. Typically, techniques using TiO bands 
to determine spot temperatures and surface areas tend to 
find larger covering fractions than Doppler imaging tech- 
niques (though see also Berdyugina (2002) J . Such differences 
in apparent spot coverage only provide a hint for unresolved 
starspots. One star for which the distribution of spot sizes is 
known in great detail is the Sun. |Bogdan et al. (19887| found 
the size distribution to be well represented by a lognormal 
function. This implies that the number of small sunspots 
is much larger than that of large spots. This supports the 
idea that there could be additional small, i.e. unresolved, 
starspots on more active stars as well. The total starspot 
coverage is of interest as a measure of stellar magnetic activ- 
ity, in order to establish the proper ratio of starspot to stellar 
plage IIRadick et al. 1998t and to obtain improved estimates 
of the total magnetic flux carried by the star (with the pos- 
sible exception of Zeeman Doppler imaging where most of 
the magnetic signal appears to come from penumbral-type 



structures, techniques of stellar magnetic field measurement 
mainly sample plage fields llSaar 19861 ISolanki "199211 1. 

Here we explore hypothetical scenarios for extrapolat- 
ing the solar spot-size distribution to activity levels typical 
for much more active stars. The basic assumption is that the 
size distribution of star spots can be described by a lognor- 
mal function, as in the case of sunspots. This assumption is 
not unreasonable since the magnetic fields on both the Sun 
and on more active cool stars are thought to be produced 
by a dynamo residing at the base of the convection zone 
l |Petrovay 2001|ISchussler fc Schmitt 2002t . From there flux 
tubes carry the field to the solar surface. The fragmenta- 
tion of these tubes during their passage through the convec- 
tion zone is thought to give rise to the observed lognormal 
distribution l |Bogdan et al. 1988^ . Lognormal distributions 
can, however, differ significantly from each other in their 
parameters. 

In order to constrain these parameters for active stars 
we investigate the possible range of behaviour between so- 
lar activity minimum and maximum and use these to ex- 
trapolate to larger levels of activity. Hence we assume that 
the processes which lead to the flux-tube size distribution 
do not change qualitatively with increasing activity. Such 
an assumption has in the past helped to reproduce, e.g., 
the high latitudes of starspots itSchiissler fc Solanki 199'2l 
ISchiissler et al. 19961 [Schrijver fc^hle 200Tt , or the pres- 
ence of active longitudes on the Sun and Sun-like stars 
[Berdyugina fc Usoskin (2003)| 
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Table 1. A selection of RS CVn stars whose spot cover- 
ages have been determined using different techniques. The first 
column gives the name of the object, the second and third 
columns the surface and spot effective temperatures according 
tojO'Ncal, Ncff & Saar (1998) The fourth column gives the in- 
clination angle, usually taken from the Doppler imaging papers 
listed in Sec. 2.1. Note that not all groups agree on the spot tem- 
peratures and inclination angles. 



Object 


Tstar 


Tgpot 


inclination 


II Peg 


4750 


3530 


60° 


EI Eri 


5600 


3700 


46° 


a Gem 


4600 


3850 


60° 


DM UMa 


4600 


3570 


55° 


HD 199178 


5350 


3800 


40° 



2 DIFFERENT MAPPING TECHNIQUES 

To our knowledge, TiO modelling has so far mainly 
been published for giant stars. One exception is 
[Saar fc Neff (1990) who calculate filling factors for 2 
As they do not give the epoch for their 



who calculate 

dwarf stars. 

observations, comparison to Doppler maps is difficult. We 
therefore limit the discussion to 5 RSCVn stars for which 
spot covering fractions have been derived from (near) 
simultaneous data using different techniques. Tab. lists 
the stars together with their effective and spot temperatures 
as well as the inclination angles of their rotation axes to 
the line of sight. Tab.|5|lists the spot covering fractions for 
the stars from Tab.Q 

The observations 
ing factors have all 
|Neff O'Neal fc Saar (1995]] 



for 



the 
been 



TiO 
taken 



fill- 
from 

lO'Neal, Saar fc Neff ( 1996)] 



and |0'Neal, Neff fc Sa ar (1998)' Due to activity cycles and 
in some cases due to incomplete phase coverage, there are 
considerable variations in the filling factors measured for 
the same star by the same group, but at different times. 
In order to be able to compare them to the spot coverage 
fractions listed in columns 6 and 8 we converted them to 
a "minimum" and a "most likely" spot coverage. This is 
described in more detail later. 

The measurements of the photometric spot coverages 
(column 6 in Tab. |2J are taken from [Henry et al. (1995)| 
[Rodono et al. (2000) | and [Padmakar fc Pandey (1999)| and 
are labelled (H-|-), R(-l-) and (PP), respectively, in col- 
umn 7. The covering fractions derived from photometry by 
[Henry et al. (1995)| are lower limits as they have used the 
maximum hght level during each individual observing pe- 
riod to represent the brightness of the unspotted star. If 
we use the brightness maxima over all of their observations, 
the areas need to be increased. While the exact increase de- 
pends amongst others on the spot geometry and the stellar 
and spot temperatures, we estimate that for a Gem total 
surface coverages of about 6% are more typical than the 
values given by [Henry et al. (1995)[ For II Peg we find that 
the surface coverage was more like 15% towards the end of 
1989 as weU as during 1992 September. [Rodono et al. (2000)[ 
give two different values for the spot coverage, one derived 
using a maximum entropy method, the other (higher) one 
with a Tikhonov regularisation. When they use the theoreti- 
cal maximum light-level inferred from TiO band calculations 
(see [Neff, O'Neal fc Saar (1995)[ for more details) , their spot 



areas increased by 15% and 20% of the total surface area for 
the maximum entropy and Tikhonov maps respectively. 

The last two columns of Table |5| give the 
surface coverage and the references for Doppler 
imaging determinations. The covering fractions 
are taken from JBerdyugina et al. (1998)1 (B98-I-), 



Washuettl, Strassmeier fc Collier-Came ron (19 98'J| (W98-|-) 
Washuettl, Strassmeier fc CoUier-Cameron (200TJ (WOH-), 
Dempsey et al (1992)[ ( D92-H)\ [Hatzes (1993)1 (H93), 
Hatzes (1995)[ (H95) and [Strassmeier et al. (1999)] (S99-h). 
Note that most Doppler maps do not give the spot coverage 
as a direct parameter, showing the stellar surface temper- 
ature rather than a spot filling factor. We were therefore 
only able to use a relatively small selection of Doppler maps 
where the authors had either given the spot coverage such 
as in as B98-I-, W98-I-/W01-I- (priv comm) and D92-I-, or 
where the maps presented allowed realistic estimates. 

Tab. 121 indicates that Doppler imaging and photometric 
light curve modelling tend to result in a smaller spot cover- 
ing fraction than TiO modelling^. Both, photometric light- 
curve modelling and Doppler imaging are prone to under- 
estimating spot areas. This is mainly because they are not 
very sensitive to rotationally invariant surface features, e.g., 
banded structures or low-level and small-scale distributed 
surface features. The current resolution limit of Doppler 
imaging is 3° to 5° in longitude, depending mainly on the 
star's rotational velocity, the spectrograph resolution and 
the signal-to-noise ratio that can be achieved. To put this 
into context, we recall that the largest sunspots have diam- 
eters of about 1°. 

In an ideal world where the spectral type of the unspot- 
ted parts of the target star were known with very high ac- 
curacy, or where the brightness and line profile of the star 
can be measured at a time when it is unspotted, the to- 
tal spot area at other epochs can be estimated. Such prior 
knowledge, unfortunately, is in general not available. Fur- 
thermore, it is unclear whether rapidly rotating stars ever 
are free of spots. 

TiO modelling is the youngest technique to de- 
termine the spot coverage and can be used for stars 
of any rotation velocity (see Ncff, O'Neal fc Saar (1995)] 
[O'Neal, Saar fc Neff (iggejl and .O'Neal, NeflF fc Saar (1998)] 
for more detail). Either the strength or the general shape of 
the TiO band heads is matched with a linear combination 
of template-star spectra at the effective stellar temperature 
and the spot temperature. In this way, the spot filling factor 
is determined. If more than one band head is observed it 
is possible to determine spot temperature and spot filling 
factor independently (provided the temperature response of 
the band heads is sufficiently different from each other) . The 
location of the spots is not recovered. The filling factors are 
weighted for the limb darkening, but depending on the loca- 



^ The surface coverage for a Gem derived by 
[Dempsey et al."{l 992) has been obtained using a variant of 
a line-bisector analysis rather than by Doppler imaging. 
^ Recent work has often combined Doppler imaging and 
photometric lightcurve modelling, showing that one and the 
same starspot distribution can reproduce both kinds of 
data. In some cases this has slightly increased the total 
spot coverage compared to Doppler maps alone (see, e.g., 
lUnruh, Collier C ameron fc Cutispoto 19951 . 
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Table 2. Spot area coverages (for the total stellar surface) of the stars listed in Table ITl obtained with different techniques. The first 
column gives the name of the object, the second column the date of the observations. Columns three to five give the coverages obtained 
with the TiO modelling technique. Column three gives the range of filling factors, while columns four and five give estimates for the 
total surface coverage derived from the filling factors. Columns 6 and 7 give the total spot coverage derived from lightcurve modelling 
and the corresponding reference. The last two columns give the surface coverage and the references for Doppler imaging determinations. 
See text for more detail. 



Object 


date 




TiO bands 




Lightcurves 


Doppler 


imaging 






ff 


coverage 


minimum 


coverage 


reference 


coverage 


reference 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


[7] 


[8] 


[9] 


II Peg 


Oct 89 


[43.. 55] 


42 


28 


13; 12/21 


H+; R+ 








Aug 92 


[36.. 50] 


37 


23 






10-15 


B98-I- 




Sep 92 


[43. .56] 


43 


28 


9; 13/23 


H+; R+ 








Jan 95 


[26.. 35] 


25 


16 


12/21 


R+ 


10-15 


B98-I- 


EI Eri 


Mar 92 


[23.. 36] 


25 


11 












Jan 95 


[<12..18] 


13 


5 






~ 3 


W98-f 




Dec 95 


[15. .15] 


13 


5 






< 7 


WOl-l- 


(T Gem 


Mar 90 


[14.. 26] 


19 


9 


2 


H+ 


2.9-5.4 


D92-I- 




Feb 91 


[27.. 33] 


28 


13 


4 


H+ 








Mar 92 


[10.. 20] 


14 


6 


3 


H+ 


7 


H93 




Jan 95 


[ 3.. 14] 


8 


3 












Dec 95 


[15.. 30] 


21 


11 


4.4 


PP 






DM UMa 


Jan 95 


[30.. 35] 


30 


15 






12 


H95 


HD 199178 


Oct 89 
May 89/90 


[16. .32] 


20 


8 






~ 6 


S99+ 



tion of the spots (e.g., a central circular spot, or spots close 
to the limb only), the area coverage can be about a factor 
of two smaller or larger than the listed filling factor. 

We have therefore calculated two estimates for the area 
coverage of the stars based on the inclination of the star and 
the possible spans of surface coverages for a given filling 
factor as shown in figure 8 of |Q'Neal, Saar fc Neff (1996)| 
The first estimate (listed in column 4 of Tab. |2J is calcu- 
lated under the assumption that the filling factor is just the 
fractional spot coverage. For the second estimate (listed in 
column 5 of Tab.|5J we try and estimate a minimum cover- 
age assuming that the spots are at disk centre where they 
produce the largest contribution. 

For both estimates we use the mean filling factor 

(//max + //min)/2 aS & Starting point (//max and //min 

are the largest and smallest filling factors observed dur- 
ing a given observing season). To obtain the minimum 
spot coverage we reduce the filling factors according to 
the graphs for disk-centre spots shown in figure 8 of 
[O'Neal, Saar fc NeflF (1996)1 This reduction is typically of 
the order of 50%, but depends on the value of the filling 
factor. 

The average and minimum covering fractions that are 
obtained in this way are in fact not the covering fractions 
with respect to the total stellar surface, as the polar region 
that is pointing away from the observer is never visible. As 
the covering fractions calculated for Doppler imaging and 
also lightcurve modelling assume that the invisible part of 
the star is devoid of structure, we multiply the mean and 
minimum coverages with a factor of (1 4- sini)/2 and so re- 
cover the values listed in columns 4 and 5 of Tab. |21 Apart 
from the conversion from filling factor to surface spot cov- 
ering, further errors in the TiO-band modelling can be in- 
troduced because of mismatches between the template star 
atmospheres and the actual stellar atmospheres. Note that 
TiO modelling assumes that stellar spectra are not affected 



by magnetic activity beyond the strengths of the molecular 
features. This is clearly a simplification and adds uncertainty 
to the spot coverage fractions deduced by this technique. 



3 SPOT DISTRIBUTIONS ON THE SUN 

Bogdan et al. (1988)| measured the size distribution of the 
sunspot umbral areas recorded at Mt Wilson between 1921 
and 1982. They found that the size distribution could be 
well fitted with a lognormal distribution of the form 

dN /diV\ / (\nA-ln(A)f\ 
dA^[lA)^'''^[ 2W^ )■ 

This is valid for umbral areas Au larger than 
Anin = 1.5 X 10"^Ai/2Q, where Ai/2@ = 2tvRq is 
the surface area of the visible solar hemisphere. Note 
that the total area of a sunspot As is the sum of 
the umbral area Au and the much larger penumbral 
area Ap. Typical ratios of penumbral to umbral area 
vary between about 3 and 5 l |Steinegger et al. 1990| 
TBrandt, Schmidt & Stci negger T990| 

|Beck fc Chapman 1993 j . In the following we assume 
that the penumbra is about 4 times larger than the umbra, 
so that As = 5Au- 

The three free parameters that have to be determined 
by observations are {dN/dA)m, i.e. the maximum value 
reached by the distribution, {A), the mean sunspot um- 
bral area and a a, a measure for the width of the lognor- 
mal distribution. If all area measurements are taken into 
account, these parameters take on the values (in units of 
lO~^Ai/20) oi {dN/d A)„, = 9.4, cta = 4.0 and (A) = 0.55. 
[Bogdan et al. (1988)| show that the same distribution can 
fit data from different cycles and that only variations in the 
value of [dN / dA)m are statistically significant. 
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Table 3. Variations in cja for a fixed {A) of 0.57. The first 
column labels the cycle activity: "all" indicates that all avail- 
able data were fitted; "max" and "min" are for data taken in 
the years bracketing solar cycle maximum and minimum respec- 
tively; "max+" and "min-|-" also include the ascending respective 
descending phase of the cycle. The second column gives the fit for 
{dN/dA)m, the third and forth columns give the values for cr^ 
and its standard deviation respectively. 



cycle 


(dN/dA)^ 




Acta 


[1] 


[2] 


[3] 


[4] 


all 


9.21 


3.95 


0.04 


max 


13.56 


4.02 


0.05 


max+ 


11.51 


4.04 


0.04 


min 


3.09 


3.75 


0.15 


min-|- 


4.99 


3.71 


0.07 



4 FITS TO THE SUNSPOT NUMBER 

DISTRIBUTION DURING DIFFERENT 
PHASES OF THE SOLAR CYCLE 

While |Bogdan et al. (1988)| show that it is possible to fit the 
sunspot umbral size spectrum with the same distribution, 
the best-fit distributions are marginally different between 
sunspot maximum and minimum. This difference between 
activity maximum and minimum is small, but a weighting 
towards larger spots could become important if we were to 
extrapolate to higher activity levels such as the activity lev- 
els observed on the stars listed in Tab.0 

In the following sections we explore what variations in 
a A and {A) are consistent with the solar data and what these 
would imply for the spot coverage of more active stars. As 
this depends on the number of free parameters of the fits, 
we separately investigate 2-degree fits with either (A) or cta 
fixed and a 3-degree fit where (A), a a and {dN/dA)m are 
all allowed to vary simultaneously. 

4.1 Fits with fixed (A) 

The initial fits involved varying only a a and (dN / dA)m- In 
the first instance we obtained a fit using the logarithmic 
form of Eq. Q Because of the logarithm, however, the er- 
rors are no longer normally distributed, so that there is no 
easy way to establish confidence limits. In order to obtain a 
reasonable estimate for the l-cr deviation, we calculated the 
value of in the region surrounding the original fit (now 
using the non-logarithmic form, i.e. Eq. 0. We hence have 
a 2-dimensional region given by {dN/dA)m and a a where 
contours of constant describe ellipses. The constant-x^ 
contours can be used to define a l-cr confidence region. For 
two degrees of freedom, the difference in between this 
region and the best fit is 2.3. The value for {dN/dA)m, <^a 
and the deviation on a a are given in Tab.|3| In all cases the 
original "logarithmic fit" lies within the l-cr contour. All fits 
listed are for {A) = 0.57. 

Note that the fits to the data taken during solar mini- 
mum when the number of sunspots on the disc is very small, 
are much less well constrained than those taken during maxi- 
mum. As a consequence, the l-cr deviation at solar minimum 
is three times larger than at solar maximum. To achieve 
statistically more meaningful fits, we also considered a com- 
bined data set of solar minimum and the descending phase of 




1 10 100 

umbral area 



Figure 1. Plot of the spot number density as a function of um- 
bral area (in units of 10~® solar hemisphere area) according to 
[Bogdan ct al. (1988) and lognormal fits to the data. The dia- 
monds show the number density at solar maximum, the triangles 
the number distribution at solar minimum. The solid lines are the 
best fit for all data sets (i.e. row 1 of Tab. |3, multiplied by the 
maximum value of the distribution at maximum and minimum. 
Also shown are the best fits to data taken during solar maximum 
(dashed line) and solar minimum (dotted line, rows 2 and 4 of 
Tab. Irrespectively). All fits are for varying a a with (A) fixed at 
0.57. 

the cycle (here labelled min-|-) and of solar maximum com- 
bined with the ascending phase of the cycle (labelled max-|-). 
While the best fits at solar minimum and maximum only de- 
viate by about 2 Acr^, the best fits at phase min+ and max+ 
differ by 4 to 5 Acr^. Some of the fits for solar minimum and 
maximum are shown in Fig. along with the spot number 
distribution at solar minimum and solar maximum taken 
from Bo gdan et al. (1988)| This shows that the distribution 
at solar minimum and maximum are indeed similar. 



4.2 Fits with fixed cja 

The fits for varying values of {A) and {dN/dA)m were ob- 
tained in a similar manner to the ones for varying a a- Here 
we fixed aA to be 4.0 and again looked for constant 
difference contours at 2.3 for the l-cr confidence limits. We 
find that (A) for the best fits varies between 0.49 and 0.58 
for solar minimum and solar maximum respectively. The ob- 
tained values of (A) are listed in Tab. 2] 

While varying aA increases the width of the lognor- 
mal distribution, varying (A) shifts the distribution towards 
larger Au- Over the range of measured sunspot umbral 
sizes both methods yield equally good fits. The implica- 
tions for small spots are, however, rather different. This 
is illustrated in Fig. |21 where the (normalised) lognormal 
distributions for fixed aA and for fixed (A) at solar min- 
imum and maximum are compared. The solid lines show 
the fits at solar minimum. The dashed line for varying (A) 
shows the shift towards higher umbral sizes, The dot- 
dashed line reveals the broader distribution obtained for a 
larger value of a a- Note that the umbral areas measured 
by [Bogdan et al. (1988)] range from about 2 to 100. In this 
range, varying a a or (A) gives an equivalent goodness-of-fit. 
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umbral area 

Figure 2. The lognormal fits to the sunspot umbral distributions 
at solar minimum (solid lines) and maximum (broken lines). The 
two uppermost lines show the fits for varying (A) , the lower lines 
for varying cr^. To bring out the small differences, we have taken 
the extreme values on the l-a contours. All fits have been nor- 
malised. For clarity the varying-cr^ fits have been offset. 



Table 4. Variations in (A) for fixed a a- The first column is as 
in the previous table and indicates the activity level. The second 
column gives {dN/dA)m and the third and forth columns give 
{A) and its l-a deviation. 



cycle 


{dN/dA)m 


(A) 


A(A> 


[1] 


[2] 


[3] 


[4] 


all 


9.40 


0.56 


0.01 


max 


13.37 


0.58 


0.01 


max+ 


11.32 


0.58 


0.01 


min 


3.45 


0.50 


0.05 


min+ 


5.67 


0.49 


0.02 



Hence the difference in the fits mainly affects the number 
of pores, dark structures that are on average smaller than 
sunspots. 

4.3 Variations in (A) and aA 

If all three parameters are allowed to vary, the l-a confi- 
dence region is inside an ellipsoid whose surface has a value 
of that is higher by 3.5 than the minimum value of x^- 
However, the data are not sufficient to constrain the three 
fit parameters very tightly and the l-a deviations increase 
by more than a factor of 5 compared to the 2-degree fits pre- 
sented in the previous section. Furthermore, the relationship 
between the (A) and aA parameters and the cycle charac- 
teristic is no longer straightforward. Going from solar min- 
imum to solar maximum in the previous sections implied 
larger values for (A) or a a- If {A) and a a are both allowed 
to vary, we tend to still get larger values for {A} (i.e., a shift 
towards larger mean spot sizes), though at the expense of 
a a that now decreases (see columns 3 and 4 of Tab.|^. This 
suggests, rather unexpectedly, a narrower distribution for 
higher levels of activity. 

When we look at the parameters within the l-a con- 
fidence regions, the picture becomes less clear, as there is 



Table 5. Variations in (A) and a a- The first column again lists 
the phase of the activity cycle. The second column gives the best- 
fit values for {dN/dA)m- Columns three and four gives the best-fit 
value of (7 A along with its deviation. The fifth and sixth columns 
list (A) along with its deviation. 



cycle 


(dN/dA) 


0"A 


AaA 




A(A> 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


all 


9.4 


4.0 


0.3 


0.56 


0.10 


max 


12.4 


3.8 


0.3 


0.64 


0.09 


max+ 


10.9 


3.9 


0.3 


0.61 


0.07 


min 


4.2 


4.6 


3.1 


0.39 


0.19 


min+ 


6.4 


4.3 


0.8 


0.43 


0.10 



overlap between the solar-maximum and solar-minimum pa- 
rameters. Most parameters suggest that the distribution at 
solar maximum is steeper as described above, but there are 
also some choices where the distribution at solar maximum is 
fiatter, but shifted towards smaller-sized spots. This makes 
it rather difficult to pick scaling parameters (see Sec. l5.3l for 
more details). 



5 EXTRAPOLATIONS TO MORE ACTIVE 
STARS 

In addition to {dN/ dA) given by Eq. the following quan- 
tities are of importance for the current analysis: A{dN/dA), 
Sn{A) and SaIa) (see also lSoianki T999ll . Sn{A) and Sa{A) 
are the integrals over (dN/dA) and A{dN/dA), respectively, 
both being normalised to their maximum values: 



Sn{A) = 


f 








r 


diV/iVtot, 




Sa{A) = 




dA' 1 






r 


A'dN/At^^, 





(2) 



Sn{A) describes the relative contribution of spots with 
area between ^min and A to the total number of spots, 
while Sa{A) gives the relative contribution of these spots 
to the total area covered by all spots on the solar or stel- 
lar surface. The latter is hence the key quantity to compare 
with stellar observations of different types such as molec- 
ular line strengths, which give a measure of the total area 
covered by spots, Atot (relative to the stellar surface area), 
and Doppler images, which provide information mainly on 
the spots above a certain size. Example plots of [dN/dA), 
A{dN/dA), Sn{A) and Sa{A) for different activity levels are 
shown in Figs |3 and |1] and will be discussed in the following 
section. 

5.1 Spot coverage for varying a a 

Having only two points to extrapolate from, we have in 
principle a large possible range of scalings. In the fol- 
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Umbral area A 

Figure 3. Plots of number and size distribution (dN/dA) (top) 
and A{dN/dA) (bottom) of sunspot umbrae as a function of spot 
umbral area Au in units of W~^Ai^2q - It assumed that rio- = 1, 
i.e., (Tyi is proportional to spot surface coverage. The distributions 
for the quiet Sun are represented by the thick lines. The next line 
towards the right shows (dN/dA) (A{dN/dA) for the lower plot) 
at solar maximum, while the lines further to the right illustrate 
the extrapolations to more active stars, up to a hypothetical star 
where about 75% of one hemisphere is covered by spots (see Tab. 151 
for a list of parameters). The distributions have been normalised 
to their values at an umbral area of j4 = 10 X lO~®Ai/20. The 
two dotted lines at ITOxIO'^Aj^jq and 3OOxlO-6Aj^y20 indi- 
cate resolution limits of Doppler imaging (see text for more de- 
tail). 



lowing, we require that aA scales with stellar activity as 
parameterised by the spot covering fraction, i.e., cta = 
(j^ -f Ao-(ylspot/^*)"" . The main open question concerns the 
choice of rio-. Plots of solar magnetic activity proxies, such 
as Ca II H & K or the 10.7 cm radio flux versus spot cover- 
age show very large scatter and it is not clear how the spot 
size distribution scales with magnetic activity. We therefore 
take the approach of choosing values of rio-, carrying out the 
analysis and from a comparison with Doppler imaging and 
TiO results deciding whether our choice is reasonable. Expo- 
nents Ua with values between 0.5 and 1.0 yield such results 
and will be discussed in the following. 

For the extrapolations shown here, (A) is kept fixed 
at 0.57, while we use the best-fit values of cta = 3.75 and 



spot area 




100000 



0.0 r '■■ i 

1 10 100 1000 10000 



Umbral area A 
spot area 

10 100 1000 10000 100000 




Umbral area A 

Figure 4. Plots of the integral functions Sjv (top) and Sa (bot- 
tom) as a function of umbral area A. The thicker lines show Sjv 
and Sa for the quiet Sun, while lines further to the right are for 
(linearly) increasing values of a a- See text for the definition of 
Sjv and Sa and Tab. 151 for a list of plot parameters. 



CTA ~ 4.0 at solar minimum and maximum respectively. The 
values for {dN/dA)m at solar minimum and maximum were 
adjusted so that the spot covering fraction at solar minimum 
and maximum were 0.03% and 0.3% respectively. Having 
preset the above parameters and picked a value for the ex- 
ponent Tier, we determine a a, the width of the lognormal 
distribution in the limit of zero spot coverage and Aa, the 
increase in the width with increasing activity. 

The solution for each activity level is then found by pre- 
setting [dN / dA)m and guessing a spot coverage ^spot/^*, 
and hence a new oa- Eq.0is then integrated to obtain a new 
spot coverage. This process is iterated until the calculated 
and input spot coverages agree. For exponents, ric, below 
about 0.7 this is straightforward with each next-higher value 
of {dN/dA)m yielding a solution with a higher spot cover- 
age. For steeper exponents, there is a threshold value for 
{dN/dA)m beyond which no solutions can be found. But be- 
low the threshold value there are generally two solutions, one 
for a low surface coverage and hence low a a, and one with a 
much higher surface coverage. This can be seen from Tab.|S] 
where for ric, = 1 the different parameters of the lognormal 
distribution are listed along with the calculated spot cov- 
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Table 6. Spot covering fractions and the proportion of spots 
above the "Doppler imaging threshold" for aA increasing lin- 
early with stellar activity (stellar activity is parametrised with 
the spot surface coverage). Rows 2 and 3 list the values for a a 
and {dN/dA)m- The forth row gives the fractional spot coverage 
of one hemisphere. The fifth and sixth rows give the fraction of 
spots that can be seen on a Doppler map with 4° and 3° resolu- 
tion respectively (note that row 5 (or 6) has to be multiplied with 
row 4 so as to obtain the fractional spot coverage that would be 
deduced from a Doppler image). 



solar: 


min 


max 












3.75 


4.00 


5.04 


10.4 


30.3 


62 


(dN/dA)m 


4.5 


38.5 


106 


106 


41.2 


20.8 




0.0003 


0.003 


0.014 


0.07 


0.28 


0.62 


^DlMtot 


0.003 


0.003 


0.010 


0.15 


0.61 


0.83 




0.006 


0.008 


0.027 


0.25 


0.73 


0.89 



ering fractions for different activity levels. The umbral-size 
distributions for the same parameters are plotted in Fig. |3 
the different lines correspond to the columns of Tab. El The 
quiet-Sun behaviour is indicated by the thickest and left- 
most line. The next line to the right is for solar maximum, 
with (TA increasing from 3.75 to 4.0 between solar minimum 
and maximum as outlined above. The two vertical dotted 
lines gives an indication of the currently achievable reso- 
lution with Doppler imaging. Only spots with areas larger 
than indicated by the dotted lines can be picked up. 

If we ask what fraction of the stellar surface is covered 
by spots larger than a given size, we have to consider the 
integral Sa- This integral is shown for the parameters given 
in Tab. |S| in the bottom plot of Fig. 21 The top graph of 
Fig- El shows the integral 5'jv. We have assumed that the 
current resolution limit for Doppler imaging is around 3° or 
4° . For circular features this is equivalent to effective surface 
areas of 340 and 600 x1O~®j4i/20- The effective area does 
not correspond to the actual area of a starspot, since the 
latter is probably composed of an umbra which produces a 
large contrast, and a penumbra, which produces a smaller 
contrast per unit area. Hence the effective area lies between 
the area of the starspot and that of its umbra. To quantify 
this effective area better we have carried out some tests us- 
ing a Doppler imaging code that is based on a spot-filling 
factor approach (see e.g. ICoUier Cameron fc Unruh 199411 . 
They show that the relative contribution of a spot umbra 
at 4500 K and a penumbra at 5400 K that is four times 
larger than the umbra is similar^. Such a penumbra-to- 
umbra area ratio is typical of sunspots (see |Solanki (2003)| 
for an overview). We therefore assume that the umbral area 
of the smallest resolvable spot is 17OxlO~®Ai/20, respec- 
tively 300x lO~®^i/20, as indicated by the dotted lines in 
Figs 13 m and |S1 This corresponds to total spot areas of 
85OxlO"''Ai/20 and 15OOxlO~'^y4i/20, i.e., spot diameters 
of approximately 5° and 6°. 

The percentage of spots that will be picked up by 
Doppler imaging is listed in the bottom two rows of Tab.|Sl 
for a A increasing linearly with the spot surface coverage, 
i.e. Ua = 1. The individual columns correspond to the lines 

This is of course just a rough estimate as the exact ratio de- 
pends on the particular line that is used for the mapping. 



Table 7. Spot covering fractions and the proportion of spots 
above the "Doppler imaging threshold" for a a increasing as the 
square-root of stellar activity (see Tab. 151 for explanations of the 
symbols). 



solar: 


min 


max 












3.75 


4.00 


4.47 


5.24 


6.6 


9.6 


(dN/dA)m 


4.5 


38.5 


155 


396 


796 


1393 


^spot / 


0.0003 


0.003 


0.015 


0.06 


0.19 


0.78 


A-Di/Atot 


0.004 


0.004 


0.005 


0.01 


0.04 


0.12 




0.007 


0.009 


0.015 


0.03 


0.08 


0.22 



drawn in Fig. 21 It turns out that for very large a a char- 
acterising wide lognormal spot distributions a substantial 
fraction of spots that are present on the stellar surface are 
indeed seen on the Doppler maps. We note that this is not 
the case for slower scaling laws {ric < 0.7), where even for 
very large covering fractions less than half the spots are 
picked up on the Doppler maps. Tab.|71lists the spot cover- 
age and pick-up fraction for Wo- = 0.5. 

Different scaling laws are also contrasted in Fig.l^where 
the spot area seen on Doppler images is plotted against the 
actual spot area. The solid lines are for Ua = 1, the dashed 
lines for ria = 0.5 and the dot-dashed lines for the interme- 
diate case of Tic = 0.75. The symbols denote data points 
from Tab. |51 where simultaneous Doppler maps and TiO 
filling factors were available. The stark difference between 
the square-root and linear scaling laws is due to the pro- 
nounced shape-change of the lognormal distribution for a 
linear increase. This results in a greater fraction of large 
spots compared to "average" spots. The number of smaller 
spots of course also increases, but this is negligible because 
of their small contribution to the total spot area. In the 
square-root scaling, the increase in the total spot area is 
achieved mainly by increasing the height of the distribution 
rather than its width as indicated by the ever increasing 
{dN/dA)m in Tab. |7| This produces a much slower change 
in the fraction of large to average-size spots so that only a 
very limited number of spots fall above the detection thresh- 
old. It is clear from Tabs |S1 and |7| and from Fig. |51 that the 
curves for ricr ~ 1 and 0.5 lie sufficiently far apart to encom- 
pass most of the data points. In fact, we find that most data 
points can be accommodated with ricr between 0.75 and 1. 



5.2 Spot coverage for varying {A) 

We now keep cta fixed and extrapolate by allowing {A) to 
vary. The assumption in this case is that as we move from 
solar minimum through solar maximum and on to more 
active stars, the mean spot size, {A), increases, thereby 
shifting the distribution towards higher values of A. The 
mean spot size increases with spot coverage according to 
(A) = {A)° + AA(A=potM,)"^. A gain, the value for ha is 
the biggest unknown. Here we illustrate the results for the 
same values as for Wo-, i.e. we assume tia = 0.5, 0.75 and 
1.0. While the fits to the solar minimum and maximum dis- 
tribution with either varying a a or {A) are very similar (see 
Fig-0) the distributions differ strongly for more active stars. 
This is illustrated in Fig.l^where the number and size distri- 
butions for increasing {A) have been plotted. Note the much 
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Figure 5. Plots of the surface coverage deduced from Doppler 
imaging as a function of actual surface coverage. The solid lines 
are for cr^ increasing linearly with the surface-spot coverage 
(rio- = 1), the dot-dashed lines are for = 0.75 and the dashed 
lines for rio- = 0.5. The lower-lying thin lines are for a resolution 
threshold of 4° , while the upper heavier lines are for a resolution 
threshold of 3° . The dotted line is plotted as a help to mark out a 
situation where all spots present on the star would also appear on 
the Doppler map. Also shown are some data points from Tab. |2] 
where simultaneous Doppler images and TiO filling factor esti- 
mates are available. These TiO estimates were used to represent 
the spot coverage of one hemisphere. The points plotted are for 
II Peg (stars), EI Eri (diamonds), a Gem (triangles), DM UMa 
(square) and HD 199178 (cross). 



Table 8. Spot covering fractions and the proportion of spots 
above the "Doppler imaging threshold" for (A) increasing lin- 
early with stellar activity, as parametrised with the spot surface 
coverage. See Tab. ISlfor a description of the symbols. 



solar: 


min 


max 












0.50 


0.58 


1.36 


3.2 


9.1 


16.7 


{dN/dAU 


5 


37 


66 


37 


15 


8 


^spot / ^* 


0.0003 


0.003 


0.03 


0.09 


0.29 


0.54 


AoiMtot 


0.001 


0.003 


0.015 


0.07 


0.27 


0.46 




0.005 


0.009 


0.042 


0.15 


0.45 


0.65 



smaller number of large spots for the "shifted" distribution 
compared to the "widened" distribution shown in Fig. |3 

The extrapolations for {A) were carried out in a simi- 
lar manner as described in the previous section, with {A) = 
{A)° + AA(^spotM*)"-*. The fixed values are a a = 4.0 for 
all fits, {A) — 0.50 at solar minimum and {A) — 0.58 at solar 
maximum. Some example parameters and results are listed 
in Tab.|H| We find that the "shift" of the lognormal distribu- 
tion towards larger A is less efficient in creating spots that 
are large enough to be picked up on Doppler maps. This can 
be seen by comparing Figs |^ and |7| where the spot covering 
fractions that would be seen on a typical Doppler map are 
plotted against the actual spot covering fraction. On both 
plots the thicker lines are for a resolution of 3°, while the 
thinner (and lower- lying) lines are for a resolution of 4°. 
Fig. Q also shows the spot covering fractions that would be 
picked up for ua — 0.75 (dot-dashed lines) and ua — 0.5 
(dashed lines). Note that the pick-up rates for ua ~ 0.5 
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Figure 6. Plots of number and size distribution (dN/dA) (top) 
and A{dN / dA) (bottom) of sunspot umbrae as a function of spot 
umbral area. The distributions for the quiet Sun are represented 
by the thick lines. The next line towards the right corresponds to 
solar maximum, while the lines further to the right illustrate the 
extrapolations to more active stars where {A) increases linearly 
with the total spot coverage. See Fig. |21 for a comparison with 
extrapolations in a a and Tab. |S]for a list of the parameters used 
for this plot. 

are almost one order of magnitude smaller than those pre- 
dicted for a linear increase in (A). Our calculations suggest 
that if, indeed, the lognormal distribution scales with mean 
spot size, then some spot clumping is needed to explain the 
relatively high pick-up rates on the Doppler maps. Starting 
from solar observations as well as from observations of ac- 
tive longitudes on rapidly rotating stars, some amount of 
spot clumping is in fact expected. 

5.3 Spot covering for varying a a and {A) 

When both, ga and {A), are allowed to vary between so- 
lar minimum and maximum, we get the somewhat curious 
situation that, for the best fits, only the value for {A) in- 
creases between minimum and maximum, while the value 
of (J A decreases (see Tab.l^J. Extrapolating from there we 
can reproduce the results found in Sects ISTTI and ^7I\ i.e., 
obtain curves similar to those plotted in Figs|3and[7| Since 
there are now two free parameters (rio- and ua) compared to 
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Figure 7. Plots of the surface coverage deduced from Doppler 
imaging as a function of actual surface coverage for scaling laws 
with (A) = {Af + AA(AspotM*)"-*. The solid, dot-dashed 
and dashed lines are for exponents, n^, of 1.0, 0.75 and 0.5 re- 
spectably. The lower-lying thin lines are for a resolution threshold 
of 4°, the upper heavier lines for a resolution threshold of 3°. The 
symbols are data points from Tab.l2lfsee caption of Fig.|3for more 
detail). 

only one in the previous sections, this is not surprising and 
we do not learn anything new. Therefore we refrain from 
discussing these extrapolations in detail. 



by Doppler imaging are actually conglomerates of smaller 
spots. The low photospheric temperature of II Peg (Table0, 
however, means that CN lines blending the analysed TiO 
band are also present in the spectrum of the immaculate 
star (Berdyugina, priv. comm.). Due to uncertainties in the 
exact temperature there are also some uncertainties in the 
above result and further such calculations for hotter stars 
would be of great interest. 

If larger spots increase in number more rapidly than 
smaller spots as stars become more active, we also ex- 
pect starspots as a whole to be more numerous relative to 
smaller magnetic flux tubes, i.e., bright magnetic elements. 
The "switch-over" between activity-bright and activity- 
dark stars seen with increasing activity level (Radick et 
al. 1987; 1990; 1998; Lockwood et al. 1992) fits weU into 
this, as does the strong increase of spot area relative to 
facular area from solar activity minimum to maximum 
( Chapma n, Uookson fc Dobias 1997| l . The fact that model 
calculations based on extrapolations from solar values do 
reproduce the switch-over at about the correct activity level 
JKnaack 199811 . supports our general approach of extrapo- 
lating the size distribution of spots from the Sun to more 
active stars. 
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6 DISCUSSION AND CONCLUSIONS 

We consider some of the consequences if spot sizes follow a 
lognormal distribution as on the sun. [Bogdan et al. (1988)| 
proposed that the passage of magnetic flux tubes from the 
dynamo to the stellar surface through the turbulent convec- 
tion zone leads to a certain fragmentation of the flux, pro- 
ducing a lognormal distribution of umbral areas (flux-tube 
cross-sections). In this picture a signiflcant fraction of the 
total starspot area of any star is in the form of starspots be- 
low the resolution limit of Doppler images. The sunspot-size 
distribution is consistent with a shift towards larger spots 
at solar activity maximum. We have used the possible range 
in the solar parameters to extrapolate to higher activity lev- 
els in different ways and have compared the resulting frac- 
tional spot areas that can be resolved by Doppler imaging. 
This has been done assuming that the starspots are ran- 
domly distributed on the stellar surface, which implies that 
most starspots resolved by Doppler imaging are single spots. 
Hence we assume that starspots do not clump. If they clump 
together as in solar active regions then Doppler images de- 
tect a larger fraction of the starspots than suggested by our 
analysis. 

Recent calculations of [Berdyugina (2002) | suggest that 
at least for the RS CVn star II Peg, spectral synthesis of 
the TiO band reproduces the observations of these molec- 
ular lines without requiring starspots in addition to those 
present on the corresponding Doppler image. If starspot dis- 
tributions on RSCVn stars are indeed log-normal and fol- 
low a sun-like pattern, this result would suggest that the 
starspots are tightly clumped and that the starspots resolved 
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